co 



Decaying Asymmetric Dark Matter Relaxes 
the AMS-Fermi Tension 

Lei Feng 1 and Zhaofeng Kang 2 '* 

1 Key Laboratory of Dark Matter and Space Astronomy, 

Purple Mountain Observatory, Chinese Academy of Sciences, Nanjing 210008, China 

2 Center for High-Energy Physics, Peking University, Beijing, 100871, P. R. China 

(Dated: April 30, 2013) 



o 

<N ■ Abstract 

u 

Q-f The first result of AMS-02 confirms the positron fraction excess observed by PAMELA, but in 

the dark matter (DM) interpretation, its softer spectrum brings a tension between AMS-02 and 
00 ■ 
CN| ■ Fermi-LAT, which reported an excess of the electron plus positron flux. In this work we point out 

that the asymmetric cosmic ray from asymmetric dark matter (ADM) decay relaxes the tension, 

Q_|. and find that at the two-body decay level a bosonic ADM around 2.4 TeV and decaying to u r + 

M-r can significantly improve the fits. Based on the R— parity-violating supersymmetry with operators 

LLE C , we propose a minimal model to realize that ADM scenario: Introducing a pair of singlets 

(X,X) and coupling them to the visible sector via LH U X, we then obtain a leptonic decaying 

^ ' ADM with TeV-scale mass. 
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I. INTRODUCTION AND MOTIVATIONS 

Dark matter (DM) is commonly accepted as a major component of our Universe in the 
present era, and its fraction in the total energy budget is precisely determined to be 26% [1]. 
Nevertheless, the most confirmative evidences for its existence all come from its gravitational 
effects, rendering its particle properties barely known. In the galaxy, DM may annihilate or 
decay into cosmic ray (CR) components, such as e^ 2 and p/p, which induce the CR excesses 
or CR anomalies. The observation of such anomalies by means of indirect DM detections 
can be regarded as a smoking-gun for the DM particle evidence. 

Recently, such smoking-guns have been triggered. In the year 2008, PAMELA reported 
that the positron fraction of the CR shows a sharp rising excess over the background, 
within a region from 10 GeV to 100 GeV [2]. This excess is extended to higher energy by 
PAEMLA and FERMI [3], in spite of a larger error bar. It is of great interest to attribute 
the excess to extra e + /e~ from DM decay or annihilation. Even more excitingly, the DM 
interpretation is well consistent with another important excess from the 2009 Fermi-LAT, 
which precisely measured the total flux of electrons and positrons and observed a flat excess 
up to the TeV region [4] [40]. On the other hand, both the PAMELA anti-proton [5] and 
Fermi-LAT diffuse gamma ray [6] data are well fitted by the pure backgrounds, which then 
stringently restrict possible dark matter explanations. In practice, the inverse Campton 
scattering (ICS) process of the injected e + /e~ always produces an associated diffuse gamma 
ray spectrum. Note that the injected e + /e~ flux is proportion to Pdm(?") 2 an d Pdm(j) in the 
annihilating and decaying DM scenarios, respectively. Thereby, to evade the exclusion from 
Fermi-LAT, the decaying DM scenario is strongly favored [7]. In supersymmetric standard 
models (SSMs), the required extremely long lifetime ~ 10 26 s and leptonic final states, can be 
realized by introducing proper R— parity violating operators [9], or dimension-six operators 
suppressed by the GUT-scale [10, 11] and further aided by some flavor symmetry [11]. 

Very recently, the long-expected AMS-02 experiment released the first result of positron 
fraction measurement [12]. It measures the fraction over a much wider energy region, from 
0.5 to 350 GeV, with an unprecedented high precision. It confirms the excess, and largely 
speaking, is consistent with the PAMELA result. Nevertheless, at the same time it brings a 
clear tension with the Fermi-LAT total flux measurement [13-16, 18]. The tension is ascribed 
to the fact that the first AMS-02 result exhibits a softer behavior than the PAMELA result at 
higher energy (As Ref. [12] emphasized, the slope of the spectrum decreases one magnitude 
of order from 10 to 250 GeV). Consequently, as we fit Fermi-LAT using a relatively heavy 
DM, the resulted positron fraction at higher energy will be too hard to fit AMS-02 (The same 
problem may also be encountered in the pulsar explanation [17, 18]). To reconcile AMS-02 
and Fermi-LAT, we may need to consider nonconventional scenarios of CR, which can give a 



relatively harder e~ spectrum than the e + spectrum. For instance, an astrophysical solution 
like hardening the primary e~ spectrum at higher energy [13, 18, 19]. 

Viewing from the particle physics, hardening the e~ spectrum can be naturally achieved 
by updating the decaying DM to the asymmetric decaying DM. In this framework, when 
the asymmetric DM decay produces asymmetric final states, say fi~T + but with vanishing 
conjugate states fi + T~ , then the asymmetric cosmic ray with harder e~ spectrum (than the 
e + spectrum) is generated. Actually, such a scenario has already been studied even at the 
time when only the PAMELA data is available [20, 21]. AMS-02 result may favor it [41]. 

This paper is organized as follows. In Section II, a model independent fit of the AMS-02 
and Fermi-LAT data is employed, assuming that the (scalar) asymmetric DM asymmetrically 
(two-body) decays to a pair of charged leptons. In Section III, we propose a minimal 
supersymmetric decaying asymmetric DM model. The Section IV is the discussion and 
conclusion. Some necessary and complementary details are given in Appendices A. 

II. ASYMMETRIC COSMIC RAY RELAXES THE AMS-FERMI TENSION 

As has been pointed out in the introduction, to well fit the high precision experiments 
AMS-02 and Fermi-LAT data simultaneously, the main difficulty lies in fitting the higher 
energy of AMS-02 (after fitting Fermi-LAT), which gives a relatively soft positron fraction 
spectrum. A way out of this difficulty is through going beyond the symmetric CR and 
assuming the injected CR is asymmetric, namely the fraction 

Generically, r + (E) has a energy dependence. Note that such a solution is thanks to the fact 
that Fermi-LAT is bind to the sign of charges. In this work, we consider the following three 
asymmetric modes from a scalar DM decay (See also [20]): 

X ->• e~/i + , e"r + , h~t + . (2) 

In Ref. [21], the fermionic DM decaying to lepton plus a charged Higgs boson or W boson 
is considered, but that may produce too much anti-protons and thus is disfavored (In fact, 
the fits are even worse than the ordinary case). 

We now employ data fitting and find out the best mode. The propagation of CR in the 
Milk Way is described by a Boltzmann equation, and in principle we can obtain the fluxes 
of CR components by solving this equation using the corresponding boundary conditions. 
However, it is difficult to analytically solve the propagation equation, due to the complicated 
distributions of sources, interstellar matter, radiation field and magnetic field. So, here we 
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FIG. 1: Top-panels: Fits for the AMS-02 positron fraction spectrum (right) and the Fermi-LAT 
e~ + e + flux spectrum (left), using an asymmetric cosmic ray originating in an asymmetric dark 
matter decaying to pT + r + . The best fit is achieved for the dark matter with mass 2.4 TeV and 
lifetime 2 x 10 26 s. Bottom-panels: A comparison to the best fit using the symmetric cosmic ray 
from DM decaying to t + t~ . 

use the GALPROP [22] package to numerically calculate the propagation of particles in 
the CR, and we modify the code to include the contribution from DM decay. During our 
fitting, we take the following parameters for the CR propagation : The diffusion coefficient 
D Q = 5.3 x 10 28 cm 2 s" 1 , the diffusion index S = 0.33, the Alfven velocity V\ = 33.5 km s" 1 
and the Halo height z^ = 4 kpc. Furthermore, the injection indexes of nucleon below and 
above the break rigidity pbr = H-5 GV are 1.88 and 2.39, respectively. Finally, as for the 
DM distribution function, we use the Einasto density profile [23]: 



p(r) = p s exp 



2 
a, 



1 



(3) 



where a s = 0.17, p s ~ 0.14 GeV cm" 3 and r s rs 15.7 kpc. 

The result is shown in Fig. 1. Among the modes listed in Eq. (2), we find that the best 
mode is X — > p~r + , with a x 2 — 90.4. The corresponding ADM mass and lifetime are 2.4 



TeV and 2 x 10 26 s, respectively. The fits for the AM S-02 and Fermi-LAT data are shown in 
the top panels of Fig. 1. As a comparison, we also display the fits in the symmetric decaying 
DM case, see the bottom panels of Fig. 1. The best fit is achieved for a similar DM but 
decaying to t + t~ , and the chi square is much worse, x 2 = 205.0. 

III. ASYMMETRIC COSMIC RAY FROM ASYMMETRIC DM DECAY 

In particle physics, it is natural to expect that the asymmetric CR originates from asym- 
metric dark matter (ADM) decay [20, 21]. However, for the sake of naturally solving the 
cosmic coincidence puzzle, namely Vt^h 2 : Odm^ 2 — 6:1, the ADM mass is expected to be 
around 10 GeV [21] rather than at the TeV scale inspired by solving the AMS-Fermi tension. 
Additionally, it is difficult to embed ADM into the simple supersymmetric standard models 
(SSMs) such as the minimal-SSM (MSSM) and its singlet (e.g., the NMSSM) or right-handed 
neutrinos extensions, owing to the robust neutralino mediated charge wash-out effect [25]. 
In this section we will first construct a simple and natural model to address these problems 
and then discuss the phenomenological aspects of the model. 

A. TeV-scale decaying ADM 

Within the conventional ADM framework [26], dark matter carries a generalized lepton or 
baryon number through proper couplings to the standard model (SM) matters. Effectively, 
this can be described by the operators such as Ot = OdmOsm(£, q) which, above some 
temperature To, establish the chemical equilibrium between the dark sector and visible 
sector, and leads to [27] 

Z>* = o, ( 4 ) 

i 

with jifc the chemical potential [27] of the particles appearing in the previous operators. In 
other words, they transfer the visible sector matter asymmetry (Assumed to be generated 
via some visible sector dynamics) to the hidden sector, which is symmetric if Ot is turned 
off. Eq. (4) indicates that the chemical potentials of two sectors should be at the same order. 
Further, for the particles in the thermal bath with temperature T, their asymmetries can 
be expressed with the corresponding chemical potentials [28] 



T 3 fj, f°° x 2 exp[- Jx 2 + {m/Tf 



n + ~ n ~ — 9~7^ / dx 

n TJo (9 + exp[-y/x 2 + (m/T) 2 } 

fb(m/T) x g b ^- (£) , (for bosons) 
ff(m/T) x p/^r- (^) , (for fermions) 



(5) 



where 9 takes 1 and -1 for a boson and fermion respectively, g denotes the internal degrees 
of freedom. The Boltzmann suppression factors fbj(m/T) indicate the threshold effects 
for heavy particles. For particles in the ultra-relativistic limit, i.e., m <C T, foj tend to 2 
and 1, respectively. Oppositely, the asymmetry of the decoupling particle can be greatly 
suppressed. Provided that only the asymmetric component leaves, we then have 

il x h 2 m x gxtxHx 

with nix the DM mass. Thereby, when the transfer ceases at a temperature far above mx, 
the conventional scenario [26], it is found that mx ~ 10 GeV is predicted to solve the cosmic 
coincidence puzzle. By contrast, when the chemical equilibrium decoupling happens during 
the period of DM entering non-relativistic, then the puzzle has to be resolved by a heavy 
DM [29, 30]. In spite of a sensitive dependence on the ratio m/T, ADM in this scenario 
still has a dynamical origin for its similar relic density to the baryonic matter's (Beyond the 
chemical equilibrium mechanism to generate DM asymmetry, there are other ways to get a 
heavy ADM [31].). 

To realize a leptonic decaying ADM, it is tempting to consider that Ot does not involve 
quark matters and moreover provides a path for DM decaying into leptons. It can be 
achieved by forcing X to develop a vacuum expectation value (VEV) vx, which breaks the 
symmetry for DM stability. However, bare in mind that here the decaying DM is extremely 
long lived, so either a very small Vx or exceedingly suppression from the operator coefficients 
is indispensable. The former, in principle is possible given some complicated dynamics, while 
the latter is inconsistent with the TeV scale ADM set up because it requires the operators to 
decouple below nix- Then we are led to the scenario where an extra source slightly violates 
the symmetry protecting DM stability. A good case in point is the R— parity-violating 
SUSY, where the tiny R— parity violations are spectator to the usual dark matter dynamics 
and only account for DM late decay. So, in the following subsection we will introduce a 
leptonic decaying ADM based on the R— parity-violating SUSY. 

B. A minimal supersymmetric model 

Now we embed a heavy decaying ADM in SUSY. Due to the neutralino mediated wash- 
out effect [25] , ADM is difficult to be accommodated in the popular supersymmetric models. 
Within the chemical equilibrium framework, Ref. [29, 30, 32] considered high dimension oper- 
ators such as Ot = X 2 LLE C with a new and relatively low cut-off scale. Here we construct a 
more economic model, even without incurring any new scales in the superpotential. On that 
purpose, a minimal model may take a form of (Actually, it resembles the supersymmetric 
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inverse seesaw model proposed in Ref. [33], where a light ADM is studied): 

W =\ l3k UL 3 E c k + (y lX L t H u X + X X SXX) 

+ (\SH u H d + |s 3 ) + W^mssm. (7) 

The corresponding soft terms are implied. Li and Lj are asymmetric, and thus Ajj& = —Xjik- 
We order i < j hereafter. The superpotential is divided into three sectors: the R— parity- 
violating sector which accounts for the ADM leptonic decay, the dark sector, of which a single 
family of lepton will be considered for simplicity (hereafter yix — > yx), and the ordinary 
NMSSM sector, i.e., terms in the second line, which dynamically generates the low energy 
scales via the singlet VEV (S) = v s . In this model the dark secfields (X,X) respectively 
carry lepton number —1 and +1, by virtue of the renormalizable transfer operator Ot = 
LH U X. However, the lepton number is explicitly violated by the first sector, so, to guarantee 
the general structure of the model, we may have to turn to other symmetries, such as 
the generalized Z% symmetry of the NMSSM. We will specifically address this problem 
somewhere else. 

We now discuss several prominent features of the model's parameter space. In the first, 
we will find that a moderately large Ax, endured by perturbativity of the model up to GUT- 
scale, is favored to annihilate away the symmetric component. So, a properly large singlet 
VEV, say v s ~ 0(5) TeV, is necessary to make Mx = Xxv s at the TeV scale. In turn, a 
large \A\ may be required in the NMSSM Higgs sector. Next, the soft spectrum is rather 
heavy because the TeV-scale ADM is also the lightest sparticle (LSP). As a consequence, 
the Higgsino and singlino masses \x = Xv s and M-g = 2kv s should be sufficiently heavy, which 
means that both A and k should take larger values. Finally, yx is an important parameter 
controlling the chemical equilibrium decoupling temperature, and it is can not be too small 
otherwise To < mx is impossible. We will discuss it more detailedly later. 

1. A Sneutrino-like ADM 

In our model, the scalar components of the chiral superfieds X and X, denoted by the 
same letters, transfer odd under the ordinary R— parity. Moreover, the dark states do not 
conserve a separated dark sector symmetry and thus the ADM candidate actually is also 
the LSP of the complete model. That is to say, ADM must be the lighter complex scalar 
from the (X,X*) mixture [42]. We now explicitly check the dark sector mass spectrum. In 
the first, the Dirac pair approximately has a mass Mx, at the TeV scale. Secondly, (X,X*) 
mix with the left-handed sneutrino z/£ after electro-weak symmetry breaking. In the basis 



[b>L,X*,X), they have a mass square matrix 



M 



( m\ y x A yx v u y x M x v u \ 
M x +m 2 x M x A Xx 

M x +m x J 



(8) 



Since the mixings with ul are suppressed by y x v u , then for a heavier m~ it is justified to 
first diagonalize the (X, X*) subsystem and get two eigenstates: 

X 1 = cos 6 X X + sin 6 X X\ X 2 = - sin 9 X X + cos 6 X X\ (9) 

with 9 X the mixing angle. The corresponding mass squares are given by 

,2 



r2 m x +m x , 



m{ 2 = M X + A 2 x T^(m x - m\f + 4(M x A Xx ) 2 . (10) 

X 1; having a lighter mass square, is the ADM (We will identify Xi\vacuum) and X\ \vacuum) 
as DM and anti-DM, respectively). It is noticed that m\ can be much smaller than M x in 
the presence of a large A\ x or negative soft mass squares. 

The mixing between X± and vl endows ADM with sneutrino properties and hence allows 
it to two-body decay alone LLE C into a pair of leptons. In addition to that, the mixing 
matters in the ADM chemical equilibrium with neutrinos discussed below. We can estimate 
the mixing angle: 

C xv L ^ 2 X - 2 (A yx cos 9 X - M x sin 9 X ) ■ (11) 

m~—m{ 

Li 

For y x ~ 0(0.1), it is can be readily arranged to be less than O(10~ A ), say by lowering sin 9 X 
and moreover taking a relatively small A yx . This helps to alleviate the need for exceedingly 
small R— parity violations. 

2. ADM at the earlier Universe 

At the earlier Universe, the dark sector enters thermal equilibrium with the plasma via its 
significant coupling to S. But such interactions conserve the dark number and thus do not 
transfer asymmetry from the visible sector to the dark sector. As mentioned before, transfer 
only proceeds after the LH U X term becomes active. Consider the neutralino-neutrino-DM 
interactions 

£ X v 3 y X U U mXil>P R Xm + C.C. (12) 



We have used H® = Yl m =i U U mXm with Xm denoting the the NMSSM five neutralinos having 
Majorana masses M Xm . Contributions proportional to Cxv L are safely neglected. Interac- 
tions in Eq. (12) lead to the following scattering 

X x v <->■ X{v, (13) 

which then maintains the chemical equilibrium between the dark states and the light species. 
The chemical potential of DM is determined to be 

Hx = -Pv (14) 

\x v expresses the lepton chemical potential and is related to the quark chemical potential. It 
is positive, so we have a negative fix and then only DM survives. Recalling that ADM should 
be about 2.5 TeV, it is reasonable to consider Tjj > 200 GeV, lying above the sphaleron 
decoupling temperature T sp . In this case, and further assume that To is above the critical 
temperature T c , we then have 

^ = J-£WT D ). (15) 

To get it we have decoupled the sparticles because we are considering a rather heavy SUSY 
spectrum. The factor kx depends on the dark sector mass spectrum, and its value varies 
between 1 and 3. We will comment on it soon later. 

On our purpose, the above equilibrium should break at a lower temperature To, when 
DM becomes semi-relativistic, namely xd = T^/mi ~ (9(0.1). As a rough estimation, To 
can be determined by equaling Txv(Td), the thermal average scattering rate of Eq. (13), 
to the Hubble expansion rate H(To) = 1.66gJ T D /M-p\- Here g* is the effective relativistic 
degrees of freedom at To-, and in our model it is at the order 100. The rate is estimated as 

T Xv (T D ) * n x (ax,v) * (o.OG^m^VV-) ^f^^ ^ W 

with gx 1 = 2 the internal degrees of freedom of the complex scalar ADM. Note that one of 
the Feynman diagrams for the scattering cross section has neutralinos in the s— channel, so 
it shows a resonant behavior as M Xm approach to m\. 

We now show some numerical results. Combined Eq. (15) with Eq. (6), it is found that 
to get a TeV-scale ADM, say mi = 2.4 TeV in our model, xd should be around 0.08, which 
is neither sensitive to kx nor the ADM mass. Moreover, despite of a strong dependence on 
mi, To does not depends on kx much. And To ~ 200 GeV is almost fixed, see Fig. 2. This 
justifies the previous assumption on T D . In addition, to decouple the chemical equilibrium for 



Eq. (13) at the proper temperature, one generically needs a sufficiently large yx- Concretely 
speaking, in our model, we need 



Ml lrr& 



\yxU um \ \ n/r2 y 2 \. 2 ~ 1(T 9 . (17) 



So, yx ~ 0.1 may be needed, even if a moderately large resonant enhancement is provided. 




0.06 0. 



FIG. 2: A plot of heavy ADM mass solution to the cosmic coincidence puzzle in our model. 

Comments are in orders. First, we have neglected other possible scattering process such 
as X\v — > tt, of which the scattering cross section is at order 0{y\ x yf). However, stops are 
much heavier and therefore there is inadequate kinematic energy (Typically E ~ mx/xo ~ 
100 GeV) available for X and u, as renders the scattering decouples much earlier. Next, 
since here the scattering is just the processes which maintains the chemical equilibrium, the 
neutralino mediated wash-out is not of concern. Finally, other dark sector states, the Dirac 
pair (X, X) as well as X 2 also store appreciable DM asymmetries given that their masses are 
near mi. The former can establish chemical equilibrium with neutrinos via Higgs mediated 
scattering. So we introduce the factor kx, having been used previously, to account for this 
uncertainty: 



V^ 9xt f{mxjT D ) 

Xi 



> — — t, — 77FT' ( 18 ) 

^ 9x 1 /(mi/To) 



with Xi running over X\p and the Dirac pair. Clearly, we have 1 ^ kx ^ 3, with the upper 
bound saturated when Xi are degenerate. 

As the Universe cools down further, the ADM annihilates away the symmetric component 
by means of DM and anti-DM collisions. We need a moderately large annihilation rate > 10 
pb [30] . For a TeV scale DM, generically this is problematic if we require \x and A sufficiently 
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small so as to maintain the petyrbativity of the model up to the GUT scale. Fortunately, the 
singlet sector is heavy (for simplicity, decoupled from the Higgs doublet sector), and thus 
the process X\X* — > H\R.2 can be resonantly enhanced by the singlet-like CP-even Higgs 
bosons H 3 in the s— channel. Here H\$ are assumed to be H u — (thus SM-) and H^— like, 
respectively. In this non-mixing limit, it is readily to calculate the cross section 



Ax (Ml\ ({\A x r 
64^ \~m{) \~^f 



°h iH2 v ~-£ -f "-— s 4- f(4mi/m 



X ) \ ttA^I ,2 



=7 . 8 x (W^y f^xy pb , (19) 

where the Breit-Wigner enhancement factor / can be tuned to larger values, e.g., 100 in 
the above estimation. Note that a large trilinear soft term A\, which is necessary to get a 
large v s , also helps to increase the above cross section at the same time. In a word, it is not 
difficult to get a sufficiently large cross section to annihilate away the ADM symmetric part. 

3. ADM today 

We now show how does the ADM leptonic decay produce the desired asymmetric CR in 
our model. Xi shares the couplings of left-handed sneutrinos vl,ii then, alone the R— parity 
violating operators LLE C , it can (only) decay to a pair of charged leptons. From the data 
fitting in Section II, the fi~r + mode is singled out. To this end, we let X123IIX1 dominate 
over other similar products (This may be arranged by a flavor symmetry) and the resulted 
operator for ADM decay is 

Aiecay ^ Cxu!^123 XifP L fi. (20) 

It is crucial to note that its Hermit conjugate part is irrelevant to the ADM decay, because 
only DM is left today. In other words, for AMD decay, the mode with charge conjugate final 
states is absent. Consequently, here we get an asymmetric CR due to the different spectrum 
between fj, and r. In addition, the ADM life time is 

When Cxv L ~ 10~ 4 or even smaller, A123 can be obviously larger than those in the previous 
works [9]. In the Appendix. A we will discuss several ways to get a very small A123. 

It is worthwhile to note that in our model the sneutrino-like ADM decay does not produce 
other significant signatures such as neutrino flux or anti-protons (But electroweak correc- 
tions on the charged particles or r hadronic decay can still induce a correlated anti-proton 
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signature [14]). The diffuse gamma ray from ICS and tau lepton decay is an exception. As 
the introduction stressed, it is stringently constrained by the Fermi-LAT gamma ray mea- 
surement. Here the bound is weakened because of the single r rather than pair r from ADM 
decay, rather than annihilation. Actually, it leaves an associated gamma ray signature for 
future observation. 

Finally we briefly discuss the prospect on the direct detection of the ADM in this model. 
The left-handed sneutrino fraction within the ADM is negligible, which implies that the 
Z— boson mediated ADM-nucleon spin-independent scattering can not be detected. Consider 
the F— term of S: \Fs\ 2 = \\xXX + \H®H® + ...\ 2 , from which we get the vertex (The Higgs 
sector is still in the decoupling limit) 

r \\ x vsm2f3sm26 x \X 1 \ 2 H 1 . (22) 

v2 

So, the SM-like Higgs boson H i = h leads to a detectable spin-independent scattering cross 
section, which can be written as a p = fp^/it [34] with n p ps m p the reduced DM-nucleon 
mass. Specified to this model, we have [35] 

m n 1 



/„ = —A Ax sin 2/3 sin 29 x 



2m,\ m^ 



-9 /'AAxA /sin 2/3 \ fsin29 x \ /2TeV\ 2 



- ° 8 x 1(r " x CiF) {"of) {-uT) {—) GeV " • < 23 > 

with m>h = 125 GeV fixed. For a p from the above parameterization to maximize the scatter- 
ing rate, the ADM has been excluded by the latest XENON100 data [36]. Of course, one can 
readily suppress this rate by turning to a much larger tan (3 and/or a smaller 9x — > (Thus 
ADM is a pure state, either X or X*). In summary, this ADM can be directly detected, 
however is difficult to be excluded. 



IV. CONCLUSIONS AND DISCUSSIONS 

Recently, the first result of AMS-02 showed the cosmic positron fraction excess. Although 
confirms the previous conclusion by PAMELA, its softer spectrum feature leads to a tension 
with the Fermi-LAT e + and e~ total flux spectrum measurement. It implies that we may 
need a somewhat unconventional extra source to resolve this tension. In this work we propose 
using asymmetric cosmic ray from asymmetric dark matter late decay as a solution. We 
find that a multi-TeV ADM (around 3 TeV) asymmetrically decays to fi~r + can indeed 
significantly improve the fit using extra source from the conventional symmetric DM decay. 
Based on the NMSSM with R— parity- violating SUSY, we construct a minimal renomalizable 
supersymmetric model to realize that scenario. The model only introduces a pair of singlets 
and involves no new scales. In addition, the ADM in our model can be directly detected. 
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In the decaying ADM framework to explain the AMS-02 and Fermi-LAT anomalies, there 
are still open questions needed to be answered. First, at the tree-body decay level, it is of 
interest to consider modes containing quarks, say X — > e~ud. After hadronization, we will 
have an asymmetric proton spectrum, and p/(p+p) < 0.5. Can it lead to an acceptable anti- 
proton spectrum by PAMELA? Next, it is not clear whether multi decay modes with proper 
branching ratios can help to further improve the fits. Last, how to distinguish asymmetric 
CR generated by the ADM mechanism from by astrophysical mechanisms, say via measuring 
the e~ spectrum, is very important. 

To end up this work, we would like to comment on the constraints of the non-annihilating 
bosonic ADM. From some astrophysical objects such as the neutron star, it is can be con- 
strained or even excluded for light bosonic ADM [38] (But exception [39]). However, heavy 
bosnic ADM has not been constrained yet [37]. 
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Appendix A: Generating small R— parity violation 

The decaying ADM life time is extremely long, ~ 10 26 s, and thus how to generate that 
small R— parity violation is challenging. In this appendix we present a simple way to explain 
its smallness. To that end, we introduce a pair of vector-like lepton doublets (L v , L v ) with 
mass around the Planck scale. They mediate lepton number violation, which is induced by 
the F— term of some singlet X (namely Fx ^ 0) to ordinary leptons. The model is given by 



c j 



! d 4 9L V L^ + J dHX^ULyFl + MyLyly. (Al) 



M P1 

Actually, the first term generates a tiny mixing between the vector-like leptons and the light 
SM leptons. The mixing is greatly suppressed by Fx/MyMp\, and as a consequence the 
induced coefficients for the R— parity-violating operators LiLjE'f, are estimated as 

^ ~ X ^MyJ^{ (A2) 

It is interesting to note that A^. as small as 10~ 20 can be achieved even if Fx/Mp\ ~ mg 
with mg ~ 1 TeV the gravitino mass, the typical mass scale of soft terms in the gravity- 
mediated SUSY-breaking. In other words, X may originate in the hidden sector breaking 
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SUSY, and even be identified with the hidden sector SUSY-breaking spurion field. 
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